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Abstract 

This work is devoted to the study of how spin texture excitations are affected by 
the presence of a static nonmagnetic impurity whenever they lie on a conical support. 
We realize a number of novelties as compared to the flat plane case. Indeed, by 
virtue of the conical shape, the interaction potential between a soliton and an impurity 
appears to become stronger as long as the cone is tighten. As a consequence, a new 
kind of solitonic excitation shows up exhibiting lower energy than in the absence of 
such impurity. In addition, we conclude that such an energy is also dependent upon 
conical aperture, getting lower values as the latter is decreased. We also discuss how an 
external magnetic field (Zeeman coupling) affects static solitonic textures, providing 
instability to their structure. 

1 Introduction and Motivation 

Two-dimensional (2D) Heisenberg-like spins models has attracted a great deal of efforts in 
the last decades. Actually, they have been applied to study a number of magnetic materials 
displaying several properties^]. For instance, the continuum version of the isotropic 2D 
Heisenberg lattice (so-called nonlinear a model- NLcrM) is useful for investigating properties 
of quasi-planar isotropic magnetic samples in the long-wavelenth and zero-temperature lim- 
its. In this case, the nonlinearity of the continuum theory supports static excitations with 
associated finite energy, the so-called Belavin-Polyakov solitons|2]. Such pseudo-particles, 
like the kinks of sine-Gordon or the 't Hooft-Polyakov monopoles of Georgi-Glashow model 

*E-mails: winder@fafeid.edu.br, winder@cbpf.br 
t E-Mail: apereira@ufv.br 



1 



have their stability guarantee! by topological features of the model rather than the equations 
of motion (for a review, see, for example Ref. 0). In fact, the presence of topologically stable 
nonlinear excitations in magnetic systems may render the latter ones interesting new phe- 
nomena. As a example, we quote the role played by vortex pair dissociation which induces 
a topological phase transition in the sample, at the temperature T BKT , even though no long 
range order is established [3]. 

On the other hand, even the purest samples contain impurities (and defects, imperfec- 
tions, etc) whose presence may give rise to important new properties. This is the case, for 
instance, of artificially dopped semiconductor materials. In magnetic materials, impurities 
may be considered for improving (magnetic impurity) or for vanishing (nonmagnetic, spin- 
less) local magnetic interaction at the positions where they were placed. However, recent 
reported results have found that the antiferromagnetic correlations around a spin vacancy 
are not destroyed, rather they appear to be increased j51 EJ Ej- Such a result has leads to 
the appearance of a new type of soliton whose energy is lower than its counterpart in the 
absence of the spinless impurity [UJ |S]. This 'more fundamental soliton' has been theoreti- 
cally [3 |Hl HH [Jill \U\ studied and observed in experiments)^. 

In the work of Ref.jHj, the NLcrM supplemented by a static nonmagnetic impurity poten- 
tial was applied to study such a new soliton. In spite of the investigation be carried out in 
the low frequency limit, its results were shown to be in excellent agreement with numerical 
and discrete lattice ones jHUZj- Actually, as presented in Ref.jHj, an attractive potential takes 
place between the soliton and the spin vacancy as long as they are sufficiently close. In addi- 
tion, such a potential is expected to give rise to oscillating solitons with definite frequencies 
around the static impurity JUj and could be taken as a trapping mechanism for solitonic 
excitations|12j. 

The discussion presented above is mainly concerned to magnetic systems defined on a flat 
plane, i.e., curvatureless manifold. Nevertheless, recent works have dealt with Heisenberg- 
like spins on curved background, e.g., cylinders, cones, spheres, thorus. In these cases, a 
number of new phenomena have been described, like the geometrical frustration on spin tex- 
tures induced by curvature and/or by non-trivial topological aspects of the space manifold, 
say, angular deficit in cones, area deficit in planes with a disk cut out, and so forth (see, for 
example, Refs. [To^ IT31 ITU Hp'] ). Actually, the study of such systems may be of considerable 
importance for practical applications, for example, in connection to soft condensed matter 
materials^] (deformable vesicles, membranes, etc), and also to artificially nanostructured 
curved objects (nanocones, nano cylinders, etc) in high storage data devices fT7\ ITS]. 

Here, we shall consider the Heisenberg spin system lying on a conical background in the 
presence of a static spinless impurity. As we shall see, by virtue of an interplay between 
geometry (cone aperture angle, a) and magnetism, soliton-like excitations experience a col- 
lective effect of impurity and geometry so that, the potential between both objects appears 
to become stronger as the cone is tighten. As a consequence, we also realize that new soli- 
tonic excitations appear exhibiting lower energy as a is decreased. We discuss on possible 
consequences of these results whenever compared to the flat plane case. Furthermore, we 
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analyse the effect of an external magnetic field on static solitons and show that such a cou- 
pling becomes these excitations unstable. 



2 The model and the influence of spin vacancy on soli- 
tonic excitations on the cone 

We shall start by considering the continuum limit of the 2D isotropic Heisenberg spins model 
supplemented by the presence of a static spin vacancy potential like below: 



H 



J 
2 



dp 2 (d t nf-{Vnf V(p) 



(1) 



with J > representing the antiferromagnetic coupling (in the static limit, d t n = 0, and 
with J < 0, the model above could describe a ferromagnetic system), while the integral is 
evaluated over a conical surface with coordinates p = (p, r) related to the usual cylindrical 



ones, r, 



by: 



r y3 

P=J, T = fr. 



In addition, n = [(sin 9 cos $; sin 9 sin $; cos 9)] is the Neel spin vector state, with independent 
variables 9(p,t) and $(p, t). Furthermore, following the work of Ref. 8\, we may write the 
non-magnetic impurity, V(p) = Vi(p), as: 
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with A = a,P/P, a is the lattice spacing parameter. 



(2) 



Then the dynamical equations which follows from Hamiltonian JIJ read like (hereafter, 
m = cos(0)): 
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A straighforward calculation also gives 1 : 
^or circular conical coordinates, (p, r), we have that: 
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VVHp) = (Pp) 1 - 1 ^ j [e p cos( 7 - \r - n\) + e r sin( 7 - |r - n|)] 5{p - pi — A) . (5) 

As may be easily checked, whenever we take the static limit in eqs. above in the absence 
of the spin vacancy (thus, V(p) = 1 everywhere in Hamiltonian (0), etc), such eqs. are solved 
by static spin textures lying on the cone with the profiles (see Ref.|16j): 

m(p) = m (p) = ^—4 , $(p) = $ (r) = r, (6) 

which represent solitonic excitations with radius R = (Ppo) 1 ^ and associated energy E = 
47r Jj3. Thus, in the (3 — > 1 limit (usual planar case), the results above recover those describ- 
ing the so-called Belavin-Polyakov solitons, presenting radius R and energy Ebp = 47r J (see 
Ref.[2 ). Then, solutions (jUJ) are the conical counterparts of the Belavin-Polyakov solitons. 

As a first result, notice that such objects present lower energy whenever lying on a cone. 
Thus, if we consider a planar surface with one or more imperfections in conical shape (cusps, 
etc) then it is expected that solitonic excitations will tend to nucleate around the apices of 
such cusps, once for presenting such a configuration a smaller amount of energy is demanded. 
This could be viewed as a geometrical pinning of solitons. 
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Figure 1: Shows solitonic profile m = mo as function of distance, r, centered at the apex cone 
(r = = p). Note how such a behavior changes as (3 decreases. Here, we have taken R = 5a. 

In addition, we should notice that the magnetization variable, m, behaves in a peculiar 
way on the cone whenever compared to the flat plane case. Figure Q shows how m = mo 
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varies with distance, r. Notice that as (3 decreases (cone is tighten) spins located from the 
apex to r = R (or, equivalently p = po), present stronger variations than those outside 
soliton radius, r > R. Physically, such an asymmetry in the variations comes about by the 
area deficit (brought about by the angular deficit): as (3 gets smaller the area on the cone, 
say, from the apex to po, also decreases, lefting minor space to spins turn from m = — 1 to 
m = 0. In the limit of very small (3, magnetization should abruptly pass from mo = — 1 
to m = 0, while the change from m = to m = 1 would be performed in a so smoothly way. 

As a first step forward, let us see how spin vacancy affects static solitons. For that, let us 
consider that the interaction of spins with such a static impurity is much less than spin-spin 
coupling. Then, we shall take: m(p) = mo(p) + mi(p) and <I>(p) = $o( r ) + ®i(p) where 
(mo, $o) are given in ©, while (mi, $i) are small corrections to the former ones, provided 
by Vj. In this case, we may linearize eqs. (@J) with respect to mi and $i for p > cr//3, 
obtaining: 



V 2 $i - 2—5 £ ■ V$i = sin T )S(p - p} 

(1 - mg) pi 



(7) 



V 2 mi + 2 



(V m ) 2 + m Vm • V 
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cos(7) 5{p - pi) . 



Now, considering that the spin vacancy is located at the soliton center, say, pi = (0, 0) 
(conical apex), then sin(7) = 0, yielding $ x = constant for eq. (JTj), while (jHJ) becomes: 



7 2 _ , o (Vm ) 2 + m Vm 



VVii + 2^ £ £ " • V mi = , 
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or still, by recalling that the 1st derivative term disappears if we change to mo = u, we get: 
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whose solutions are the trivial one, mi 



(1 -u 2 
and: 
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mi(m ) = ci(l - ml) + c 2 



m o + -(1 



mo) In 



1 + m 
1 - m 



with ci and c 2 being constants to be determined by the solitonic characteristics. Indeed, by 
demanding that as p — > oo, m = m + mi ->lwe must impose c 2 = 0. Then, the complete 
solutions in the presence of the spin vacancy are the following: 



mp = m , $p = $c 
mi = m + ci(l — wig 



$ 7 = $ c 



(10) 
(11) 



Actually, we should notice that the presence of the nonmagnetic impurity support both of 
the two solitonic excitations above, similarly to what happens even in the usual plane case, 
Ref.jS]. Their energy read as below: 



J r 2 ^ 
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(12) 



where j stands for soliton type, say, j = I, P, while the lower cut-off in the p-integral, a 13 / f3, 
takes into account the effect of Vj(p)-potential. Evaluating these integrals over the cone, we 
finally obtain (recall E = Att Jf3): 



Et = En 



Ep = En 



+ (Ra) 2 ^(l - Cl ) 
{RW + a 2 ?) 2 
R 2 ? 



(R 2 ? + a 



2f3\ 



(13) 
(14) 



Now, minimizing Ej with respect to soliton size, R(a), we see that < C\ < 1/2, analogously 
to its counterpart found in Ref. [Bj. 

How such energies behave as function of soliton radius, R, is displayed in Figure El 
Similarly to the results presented in Ref.|H], /-type appear to be more fundamental than 
P-solitons. Indeed, for (3 = 1, our results exactly recover their ones. Furthermore, it is 
worthy noticing that as (3 gets smaller the associated energies to both kinds of excitations 
also decrease in such a way that the gap between Ej and Ep turns out to be larger (see 
Figure 121). Since the observation of these solitonic excitations were recently reported we 
claim that similar experiments concerning spins textures on conical support could determine 
the energy gap increasing predicted here. 

Actually, in the preceding analysis we have considered that the nonmagnetic impurity 
is located at the soliton center, say, p = (0,0) (the apex of the cone). However, following 
the work presented in Ref.|H], we may wonder whether the effective potential between them 
would behave as long as they were placed one apart other. The way to perform that is 
analogous to that presented in Ref. jH] , and we shall not repeat them here (we refer reader to 
such a reference for further details). After some calculation, we get (the soliton is centered 
at the apex, p = (0, 0) = r): 



U, 



eff 



R^aHri - a?) 



\l{4-aP) 2 + rf[{4 - - R^] yVo + a?) 2 + r 2/3 [(r£ + a/ 3 ) 4 - R^] 

1 / (Raf \ 



+ 



4 \{rf + 



(15) 



where pi = r^/ (3 labels the impurity position. Notice that the singular points of U e s appear 
at (ro)J = (R 13 =F a 13 ) 1 ^. For example, if (3 = 1 and R = 5a, we then have (ro)g = 4a and 



6 



1 2 



1 O 



2- 




toeta=0.3 



1 O 



Figure 2: Displays the energies (in units of J-parameter) associated to I and P solitons, Ei and 
Ep, represented by dashed and solid lines respectively. Note that as the cone is tighten, both 
energies appear to decrease in a such a way that the gap between Ei and Ep becomes greater, 
keeping always the more fundamental characteristic of /-type excitations. Here, we have taken 
a = 1 and c\ = 0.1. 



( r o)s = ® a - m Figure © we plot U e s versus r for some values of (3. In general, the potential 
is attractive for r < (R 13 — a /3 ) 1 / /3 while for separation r > {R? + a?) 1 " 3 its appears to be 
repulsive. Such a behavior was already obtained in the work of Ref. jHj, for the planar, say, 
[3=1 case. 



What our results bring as novelty is the fact that as long as the cone is tighten (/3 gets 
lower) then the range of the attractive potential gets shorter. As a consequence, if we fix 
R = 5a, and extrapolate our results to the discrete lattice (with spacing a), we conclude 
that for f3 < P cr (P C r = hi(2)/ ln(5) ~ 0.43, in this case), we get r < a what implies that or 
the soliton is pinned on the impurity or they are quite apart one from another thanks to the 
repulsive potential. On the other hand, keeping the R- value above, whenever /3 cr < /3 < 1, 
we have the possibility of oscillating solitons around the spin-vacancy. Actually, let us take 
U e ff above, and consider small displacements impurity-soliton (r ). Expanding it like below: 



where we have taken into account the type of the excitation, according to (jl3H14|) . Now, 
looking only for the harmonic potential (say, the quadractic term) it is easy to show (we 
refer the reader to Ref.|TDl for details) that the solitons oscillate around the spin vacancy 
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Figure 3: Exhibits the behavior of U c r (in units of J-parameter) as function of the separation 
soliton- vacancy, r*o (above R = 8a and a = 1). Note that for short distances the potential is 
attractive while for large separations it turns out to be repulsive. Notice also that such an attractive 
interval decreases as (3 gets lower values, indicating a stronger confinement (pinning) of the soliton 
whenever /3 is sufficiently small. From left to right, we have j3 = 1, 0.7, and 0.4, respectively. 



with frequencies given by: 



,2/3 



+ {R 2P + af) 



,2/3 
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, v / 5[i? 4 ' 9 - (2a) 4 / 3 ] i(R?P + a 2 / 3 ) 2 
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(17) 



where a,j = ai,ap (a/ = 1.01a and ctp = 0.23a, see Ref.jE], for details) while c = 2JSa/h = 
Ej/Mj, where Mj is the soliton rest mass (for details, see Ref.|19j). 

How ujj behaves as function of the type (J or P) and soliton size, R, for some values of 
/3-parameter is shown in Figure (jlj). As expected, and earlier reported in Ref . [TU] . as soliton 
becomes larger, their frequency modes decrease. In addition, J-solitons oscillate faster than 
P-type ones. The novelty brought about here is that, as conical support is tighten the exci- 
tations appear to oscillate faster. In addition, notice also that for a given /3 the soliton size 
should be R > 2 1 ^a (in the discrete lattice). Thus, we see that for (3 — 1 we should have 
R > 3a while for (3 = 1/2 such a radius should take values > 5a. However, despite such a size 
increasing, the geometry of the support deeply affects frequencies, as displayed in Figure (jlj). 



3 The Zeeman coupling and the instability of static 
solitonic solutions 

In this section, we shall briefly discuss how the coupling of spins to an external magnetic field 
yields unstable static spin textures. Introducing the term (figB ■ n)V(p) = « gBomV(p) 
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Figure 4: Shows solitonic frequencies, uij (in units of c/a), around the static impurity as function 
of soliton radius, R (here o = 1). Notice that /-excitations (dashed lines) present higher frequencies 
than P ones (solid lines). Furthermore, note that as long as the cone is tighten, these excitations 
oscillate faster (the two lower curves are concerned to f3 = 1 while the upper ones are related to 
= 0.5 ). 

into Hamiltonian (Q), where B = B z (B is homogeneous), the counterparts of dynamical 
equations ()3I4|) read like follows: 
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where ^1 = J/gfiB is the magnetic length. 



As we have done in the preceding case, supposing that J >> g/i B (spin-spin coupling is 
much stronger than Zeeman interaction) then, we have that eq. ((7|) is still valid here, while 
(jSJ) is modified by the magnetic field. With the same assumptions taken in the preceding 
analysis (impurity at the soliton center, cylindrically symmetric solutions, etc), we conclude 
that $i = constant remains the solution for its respective differential equation also here. 
However, eq. now gets the form (using mi = mi (mo) = mi{u) analogously to eq. 0): 



d 2 mi 2 mi 



du 2 
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(20) 
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Contrary to eq. (JHJ), the above one does not present analytical closed solutions for ar- 
bitrary f3, as far as we have tried. Indeed, by using MAPLE8 we have obatined closed 
expressions for some (3 = l/N- values (N positive integer). Namely, we have noted that such 
expressions becomes more complicated as N is rised. By virtue of their length, we shall not 
write them here. Rather, in Figure (j5J) we plot m = + mi as function of distance, r. 
The main point to be noticed here is that no good mi-functions solves eq. (J2(Jj) . what may 
be realized by noting that m blows up as r — > oo, yielding divergent associated energy. Ac- 
tually, in the very beginning we have defined m = cos(6 l ), but the results plotted in Figure 
(0) explicitly shows us that for r > L c (depending upon (3) m- variable takes values > 1, 
what is a contradiction. Thus, we conclude that whenever Zeeman corrections are taken 
on m solitonic profile these excitations appear to get unstable and collapse. Furthermore, 
once that such a blowing up character of m-profile comes about by turning on 5-field and 
appear to become stronger as long as the cone aperture gets smaller, we may realize that 
by decreasing /9-parameter we can increase the net magnetic effect on solitonic excitations. 
Such a result can be viewed as a geometrical accentuation of the magnetic field experienced 
by solitons. 




Figure 5: Displays m solitonic profile versus distance, r. Without external magnetic field, m = mo 
and m — * 1 as r —* oo. However, whenever B is turned on, such a profile present values > 1 
for r > L c , what contradicts m definition and implies that static solitons becomes unstable. In 
addition, note that such a critical size, L c gets smaller values as (3 is lowered. 

Notice, however, that we cannot state that they shrink to point-like objects, since the 
impurity introduces a characteristic length and no smaller size than spin vacancy is possi- 
ble (recall that the impurity was placed at soliton center). In the discrete lattice case, we 
would expect that Zeeman interaction should force solitons to shrink to the smallest size: 
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the lattice spacing a. Notice, however, that this is pointed out here as a possibility (following 
the conclusions presented in Refs. |13[ IT3]) once our analysis is strictly valid only for r > a. 
Indeed, a similar problem was earlier considered in the work of Ref . |14j . There, the authors 
explicitly state the impossibility of carrying out their analytical analysis for the highly com- 
plicated diferential equations. 

Even though we could not provide here a rigorous proof of whether the solitons nec- 
essarily decreases their radius to spin vacancy size in order to compensate the instability 
brought about by Zeeman interaction, we should mention that the present problem is im- 
portant in connection with dynamical solitons, for instance, in the lines of Ref. . Thus, 
considering precessional oscillatory modes of a soliton around B, their associated frequencies 
could be determined and how geometrical features affect such modes. Such an issue is under 
investigation and results will be communicated elsewhere |21j. 

4 Conclusions and Prospects 

We have considered the problem of classical Heisenberg spins on a circular conical geometry 
subject to both the potential of a static non-magnetic impurity and an external magnetic 
field. Whenever spins weekly interact with the impurity alone, a number of results were ob- 
tained and compared to their planar counterparts. In general, our results recover planar ones 
as long as /3 = 1, and for arbitrary cone aperture they describe how solitonic spin textures 
are affected by a static spin vacancy in such a geometry. Among other, we have seen that the 
impurity presence now allows the appearance of two kinds of solitons: I and P-type ones. 
In this line, one interesting thing realized in conical geometry is that the energies of both 
types get smaller values whenever the cone is tighten while the energy gap between them is 
increased. Since both kinds of excitations have been recently observed in experiments, we 
could think of similar experiments dealing with Heisenberg-like spins on conical supports in 
order to determine the results predicted here. 

Since NLcr model has been also applied to study weekly self-interacting two-dimensional 
electron gas (2DEG) [22], our present study may also be relevant, for instance, in connection 
with Quantum Hall Effect in curved, say, conical, geometry. 

An analysis of how external magnetic field influences the structure of dynamical solitons 
is in order. For instance, how their internal precessional modes and sizes are influenced by 
this field are some of the issues we have been addressing. Results concerning such a study 
will be communicated elsewhere |2Tj . 

As prospects for future investigation, we may quote the introduction of more impurities 
in the system, for example, as presented in the works of Refs. |lH |2H]. There, analytical 
and numerical/simulational analysis of the model may indicate how collective spin vacancies 
effects are read by solitonic excitations. 

Concerned with the current topic of magnetic nanostructured objects, how nanomagnetic 
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cones [T7| HH] (among other curved structures), with and without vacancies and magnetic 
field, affects the magnetization and hysteresis loops are some interesting points to be studied. 



Finally, in connection with Field Theory/High Energy Physics and related topics, a 
problem that could lead to interesting results is that of a charged particle interacting with a 
magnetic monopole in a conical background. Perhaps, in this framework, bound states be- 
tween them could be observed, scenario quite distinct from its flat space counterpart, which 
does not present closed orbit configurations between the m [23] . 
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